SOME SOLUTIONS TO EQUATIONS OF 
MOTION IN THE EQUATORIAL 
REGIONS 
R. M. JONSON 









Cm— I 

S3 

V/v 

' :5 i! 

1 i -i... . . • -• : 






I ii>rar\ 

U. S. Naval Postgraduate School 
Mont ore v, (California 



SOLUTIONS TO B^SiaU'S OF KOTlCti 
III TH3 3 JJATCRLuL 33GI0KS 



S. M. Jonson 



S0I-C3 SOLUTIONS TO EQUATIONS OF MOT I Oil 
117 THE E-^JuTGRIAL REGIONS 



t>y 

Bussell Martin ^Tonson 
Lieutenant, United States Lavy 



Submitted in partial fulfillment 
of the requirements 
for the degree of 
MASTER OF SCIENCE 
III AEROLOGY 



United States Naval Postgraduate School 
Monterey, California 

1952 



This work is accented as fulfilling 
the thesis requirements for the degree of 

KAST3R 0? SGIKjGS 
III iiSROLOG-Y 



from the 

United States naval Postgraduate School 



PREFACE 



This investigation v/as conducted at the United States ..aval 
Postgraduate School, Konterey, California as the thesis requirements 
for the degree of Master of Science in aerology. 

For help and advice received in its preparation the author is 
indebted to Associate Professor Gr. J. Haltiner of the U. S. naval 
Postgraduate School. 



do 



TABLE OF CONTENTS 



Page 

CERTIFICATE OF APPROVAL i 

PREFACE ii 

TABLE OF CONTENTS iii 

LIST OF ILLUSTRATIONS iv 

TABLE OF SYMBOLS ALT) A33REV I AT I Oil S v 

CHAPTER 

I. INTRODUCTION 1 

II. MOTION AS A FUNCTION OF y ALONE 4 

III. SPATIALLY INDEPENDENT MOTION WITH VARIABLE 

CORIOLIS PARAMETER 11 

IV. SPATIALLY INDEPENDENT MOTION FOR A CERTAIN 

PRESSURE DISTRIBUTION 13 

V. CONCLUSIONS 1 6 

BIBLIOGRAPHY 17 



(iii) 



LIST OF ILLUSTRATIONS 



Page 

Figure 1. Sea level streamlines for Uq = 5 nips, 

v 0 = 5 mps and for uo = 5 mps, v 0 = 0 at y = 0 6 

Figure 2. Sea level streamlines for 

u 0 * -5 mps, v Q = 0 at y = 0 7 

Figure 3- Sea level streamlines for 

Uq = -5 mps, v 0 - 0 at y = -500 km 8 

Figure 4. Sea level streamlines for 

u 0 - -5 mps, v 0 = -2 mps at y = 500 km 9 

Figure 5* Sea level streamlines for 

u 0 = -5 mps, v 0 - -5 mps at y = 0 10 



(iv) 



I 



table of symbols a.d ^lbbsviatiohs 

> Coriolis parameter » 2 u> sin 

f Density of air 

t Time 

cj> Latitude 

OJ Angular speed of the earth = 7*29 x lO - ^ sec - ^- 
x Coordinate axis - Last 

y Coordinate axis - Horth 

u Easterly velocity component 

v northerly velocity component 

a Badius of the earth = 6378 km 
T Period of pressure variation = 12 hours 

Ci.Cg Constant of integration 



(v) 



I. INTRODUCTION 



In mid-latitudes the geostrophic and gradient solutions of the 
equations of notion generally give a good approximation to the true 
wind. However, in the equatorial regions the acceleration terms, both 
centripetal and tangential, are frequently of the same order ox magni- 
tude as the pressure and Coriolis forces; hence the gradient and 
geostrophic solutions are not generally applicable. 

To render the equations integrable various simplifications have 
been used, Potable among these are the assumptions that the flow is 
(a) horizontal, (b) frictionless, (c) homogeneous, (d) non-divergent 
and (e) steady* Some investigators have neglected the Coriolis force, 
while others have either treated it as a constant or replaced sin i> 
by p • Practically all authors assume (a) through (c) above as v/ill 
be done by all discussions in this paper* 

In addition to the above assumptions, Grimes [$] let motion be 
independent of longitude* He allowed for variation of the Coriolis 
force by assuming sin <p = j> with as a function of y* His solutions 
gave particular streamlines which were determined from initial values 
of velocity and vorticity. Pressure distributions were also determined 
from these solutions; however, it is generally preferable to derive 
motions from known pressure distributions and initial velocities* 

Grimes obtained isobaric patterns which were similar to mean maps in 
the Indian and South Pacific oceans but pressure gradients were less 
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than those usually observed. Crossley [ 2 ] found this solution to 
apply only to mean motions and not to the instantaneous flow pattern. 

Orossley [ij removed the restriction on independence of longitude 
and corrected the fault that all isobars crossed the equator at right 
angles in the above solution. The same author [3j obtained a solution 
in spherical coordinates. 

Crimes 1 original assumptions less the requirement of non-divergence 
were used in Chapter I of this paper. 

Forsdyke fhj took a constant Coriolis force with motion at any 
instant everywhere the same, that is, the space derivatives do not 
appear in acceleration terms. He assumed a pressure distribution of 
the form p ® p 0 +A(t)y for various functions A(t). 

By expressing the resulting motion as geostrophic plus ageostrophic 
components, he suggests a method for relating wind to the pressure dis- 
tribution in the tropics. 

In Chapter II an attempt is made to generalize Porsdyke’s solution 
by making less restrictions on the Coriolis parameter and assuming a 

pressure distribution of the form p - p 0 + x sin 2H t . 

T 

In Chapter III a pressure field of the form p r p 0 +Ae~ e ^ 3(l-e“ t )y 
where A, B, € and Y are constants is discussed. This pressure field was 
suggested hut not developed by Forsdyke (VJ • 

More recently, Schmidt [?J concentrated on the kinematics of flow 
hy starting with Rossby's vorticity equation and assuming the dependent 
variables to be functions of y only. He obtained values for u in terms 
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of y assuming v known and plotted streamlines for £ monsoon current, 
high pressure at the equator (doldrums), and a line of convergence 
stretching in an east-west direction along the equator. Some of the 
streamlines obtained in Chapter II of this paper closely resemble 
Schmidt ' s. 

G-rimes [6j using his original assumptions less independence with 
longitude, obtains a set of pseudo-geostrophic equations where stream- 
lines are absolute vorticity trajectories. They are parallel to 
isobars of "dynamic pressure" (P - , where - u‘ + v^), 

The difficulty in this approach is the necessity for drawing accurate 
isobars of dynamic pressure, v/hich must be drawn at intervals of one- 
fifth of a millibar. This is particularly difficult where reduction 
to sea level is necessary. 






( 3 ) 



II. i.OTIGi: ..3 x, FJ:;CTIu» OF y ixLOi'iS 



The equations of horizontal frictionless notion are: 



^ - 2 aJv sin <P - - — 
dt f* 3.x 



( 1 ) 



— f2a?u sin <t> 


_ _ 1 2n 


dt 





In this solution steady-state is assumed with motion independent 

of longitude « iJL = 0). pressure is a linear function of latitude 

3x d x 

(y) alone v/ith isobars equally spaced and parallel to the equator. 

The Coriolis force is treated in a similar manner to Grimes £5j hy 
letting sin <f> m $ - y/a. Thus the Coriolis parameter is Ky where 
K - 2a2/a. 

'with these assumptions the equations of motion become: 



du 
v — 

BJ 



iiyv 



v 

By a/ 



Equation (3) readily integrates to 



u * i Ky 2 + u 0 



(3) 

(*0 

(5) 



Substitution of this value of u into equation (4) and integrating 
yields the following expression for v. 
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( 6 ) 



T 2 *v 0 '+2c( (p-Po) - i I2 2 - Kuo y 2 

u 0 , v o» Po being values of velocity and pressure at y - 0* 

numerical values of the slope of the streamlines were found at 
every $00 kilometers with varying initial conditions and a constant 
pressure gradient of 1.5 millibars per $00 km and a specific volume 
of 860 cm3/gm. 

These streamline patterns (see Figures 1-5) resemble observed 
mean motions in the Indian and Southern Pacific oceans. Values of 
v for y > 1500 km v/ere imaginary indicating steady-state is no longer 
valid at that distance from the equator with the given initial con- 
ditions. 

Most of the faults of G-rimes* and Crossley's solutions still 
remain in this simple pattern, the most important being that steady 
state is not generally applicable and therefore solutions apply only 
to mean motions. It should also be noted that the determination of 
the streamline patterns from even these simple pressure distributions 
is a laborious process. 

However, one improvement is the removal of the restriction of no 
horizontal divergence. Zones of divergence and convergence are known 
by observations to exist in the eopiatorial regions. 
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Sea level streamlines for u 0 = 5 mps, 
v o = 5 mps and for Uo = 5 mps, v 0 = o at y - 0 

Pig. 1 
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Fig. 2 
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Sea level streamlines for 
-5 mps, v 0 * 0 at y = -500 



Fig. 3 
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Sea level streamlines for 
u 0 = -5 zips, v 0 ” -2 mps at y * 500 km 

Fig. 4 
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Fig. 5 
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III. SP^IaLLY IPJSPaD-iTT IIOTIOL .,ITH 
VARIABLE CORIOLIS PiJRAI^IETSR 

In this approach Porsdyke^ [4j assumptions were used hut the 
Coriolis parameter was varied to secure a more general solution. 
Motions were considered uniform which is probably valid over areas 
five degrees square. The pressure attribution p = Po+PlX sin ~~ 
was assumed. 

The equations of motion with the above pressure distribution are: 



— sin 2d Lk ( 7 ) 

e f * 



( 8 ) 

dt 



Differentiating ( 8 ) with respect to t, we obtain 



d 2 v . 

d7 ' 




dA 

dt 



moreover 



A = 



sin 



4 



-7- = 2 ^cos <f> 
dt dt 



0 = y/a; and 



(9) 



d^ m _1 dy/dt * v/a. Since cos^ varies only 20$ from 0° to 20°, we 
dt a 

assume cos <fr - 1 and it follows that. 



dA - 2 u?v (10) 

dt a 

Substituting the expressions for du/dt, u, and d^/dt from equations 
( 7 ). ( 8 ), and ( 10 ) respectively, we get 
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jl _ 4s g _>>f» ^ 

its* «7) <tt ^ T 

To solve this differential equation a mean value of v(v) was 
assumed giving a second order equation with constant coefficients. 

The solution with constants of integration evaluated for the 
"boundary condition, v = 0 at t = 0 and Cj = 1 is 



a»v + -£ 

v/s e ** t(-l- 8 ) & 

■+ A Sin + 3 Cas lz± 

r ~r 



>j) U '■Itv'y 7 ' Cl))\ £ 
cl A 



a = 



\ UJV £ 

_fu> V i- \fuj*V - <C^ 1/3 /ZA 

V oTa 1 * / c 

/ > : — TZ.) 1 -/ 

[ an^-v^v - Bz 2 l I e, 

v a). 3 - ' 

- ^ -f -5 -£-!?■ J/n «?JL? 

A ~r ~T A ^ T 



where 









and 



5 = 



d ijz* t — r^* 



The special case for which v = 0 gives Forsdyke's solution. Assigning 
normal numerical values led to imaginary values of u and v and hence 
the analysis was carried no further. Some numerical constants will 
lead to real values for u and v; however, the results do not conform 
to observed patterns. 
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IV. SPATIALLY I-.D3?2L!D3l. r T K0TI01I ?OR A 
CERTAIN PRESSURE DISTRIBUTION 



Forsdyke OJ suggested a pressure distribution of the form 
p > p 0 i- Ae~^ x Etl-e^^Jy - where A, B, 6 ,f are constants t he 



applied to the spatially independent constant Coriolis parameter 
solution. He states this solution might he applied to curved isobars. 

This statement appears to he inconsistent for at any time d, the 
pressure is a linear function of x and y giving an isoharic pattern 
consisting of straight parallel lines with their orientation changing 
with time, at t = 0 the isobars are a function of x and perpendicular 
to the equator, ns t increases and approaches infinity the isobars 
rotate and become parallel to the equator and are functions of y alone. 

The constants A and B are the pressure gradients at times t = 0 
and t - 00 respectively. The constants 6 and t determine how rapidly 
the pressure gradients change in the x and y directions. 

In applying this pressure distribution fc and \vere made equal 
for simplification. 

The equations of notion with Jorsdyke's assumption and the above 
pressure distribution become 
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2s. = ->u - 2(i- e -"* b 

f 



(13) 
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The general solution of this set of differential equations 



consists of the complementary functions 

U. - strict — Cl £05 
V - C, s/n T\~6 + Cos ?\ -t 



and the particular integrals 



_ JL - {&l± Ax) pi* 
~ e7\ 






AT) -BX p -V£ 

e(r^J 



If and ^2 are evaluated at the boundary conditions u = u Q , v = v Q 
at t - 0 and the values of geostrophic wind 



Uq- - J_ B Q _ 

*/ » _L df~JLp- z ^ 

* ^ dX ^ e 






are inserted we have 



, retsM+d s/0 )n y L- §2 
u LfO>v> J L p ptf+TtX 



CCS 

(14) 



-f 



<43 UV>' 



- ?2l v? 

»V>vV ' w W3 * 



v-'H- caS/H 



(n L 

W, 





(15) 
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If a - 0 and. yf - oo the pressure distribution becomes 
p ■ p 0 i- 3y, which is the same as used in Chapter II of this paper 
and also in one of Porsdyke's [b] solutions for a particula.r time 
variation of pressure. 

Substituting A. = 0 and t - 00 in (14) and (15) we get 



u = v 0 sin A t (uq - ^A.) cos \ t 






(16) 



= (— u 0 ) sin ^ t + v 0 cos A t 



(17) 



Equations (l6) and (17) are the sane as Forsdyke obtained when 
the substitutions u^ = u 0 - Ug and vl = v 0 - Vg are made since 
Ug ® and vg * 0 for A = o , ^ ~ 00 . 

Equations (lo) and (1?) differ 'from (3) and (4) in spite of the 
same pressure distribution, This is consistent since different 
restrictions were placed on the motion; namely, constant Coriolis 
parameter and independence of space for (l6) and (17), while steady 
state, independence of longitude, and a variable Coriolis parameter 
in (3) and (4). 

numerical evaluation of equations (lo) and (17) coula be done by 
assigning values to the constants A, 3, and y consistent with observed 
pressure patterns and then computing u and v as a function of latitudes 
and time. This evaluation was not undertaken here for lack of time. 
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7. CONCLUSIONS 



The equations of motion nay be solved for equatorial regions 
only with rather restrictive assumptions even for relatively simple 
pressure fields. 

Although the solutions in this paper give velocity fields 
resembling observed mean motions, and one commonly used restriction 
has been lifted, namely, the requirement non-divergence, they appear 
to be of little practical value in routine weather analysis. In 
addition to other limitations the necessary computations are 
laborious and time consuming. 
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